Focused stimulation in a medical stimulation device

ABSTRACT

A medical stimulation device such as a cochlear implant configured to provide stimulation of one or more spatially-restricted contiguous portion(s) of the spiral array of auditory nerve fibers in the cochlear (“discrete stimulation regions”). Each discrete stimulation region is defined by the constructive and/or destructive interference of stimulating and limiting signals simultaneously applied to electrode channels of an implanted electrode array, the stimulating and limiting signals being determined based upon transimpedance measurements of intracochlear electrode channels of the implanted electrode array representing specific spread functions of an individual recipient. The stimulating signal is preferably applied through a targeted electrode channel; that is, one or more successive electrodes which is/are adjacent to the discrete stimulation region. The targeted electrode channel is selected to represent a sound based on the outputs of a sound processor to stimulate neural activity in the discrete stimulation region to thereby cause a percept of the represented sound. The size of the discrete stimulation region is defined by the limiting signal(s) applied to electrode channel(s) other than the targeted electrode channel, and which negate(s) current spread which would otherwise occur in response to the stimulating signal.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims priority from U.S. Provisional Patent Application No. 60/675,860 entitled, “Phased Array Stimulation,” filed Apr. 29, 2005, which is hereby incorporated by reference herein.

BACKGROUND

1. Field of the Invention

The present invention relates generally to medical stimulation devices, and more particularly, to focused stimulation in a medical stimulation device.

2. Related Art

There are several types of medical devices that use electrical signals to activate nerve, muscle or other tissue fibers in a recipient (also referred to herein as a user, listener, patient, etc.; “recipient” herein) to stimulate an activity. Such medical devices are generally referred to herein as medical stimulation devices. Usually the activity is invoked to compensate for a deficiency in the recipient. For example, stimulating prosthetic hearing devices have been developed to compensate for loss of hearing.

Several types of prosthetic hearing devices provide electrical stimulation to aid recipients who have a hearing deficiency. For example, cochlear™ implants (also referred to as cochlear™ devices, cochlear™ prosthetic devices, cochlear™ implants, and the like; simply “cochlear implants” herein) apply one or more stimulating signals to the cochlea of a recipient to stimulate hearing. Another example is an auditory brain implant that delivers electrical stimulation to the auditory brainstem nuclei of a recipient to stimulate hearing. For ease of description, the present invention is presented in the context of a prosthetic hearing device, namely, a cochlear implant. However, it should be appreciated that unless stated otherwise, the present invention is applicable to any medical stimulation device now or later developed.

Cochlear implants typically include a sound transducer that receives incoming sound, and a sound processor that converts selected portions of the incoming sound into corresponding stimulating signals based on an implemented sound encoding strategy. The sound processor transmits the stimulating signals along an electrode array implanted within or adjacent to the cochlea of the recipient.

Cochlear implants exploit the tonotopic organization of the cochlea by mapping audio energy in specific frequency bands to deliver stimulation at corresponding locations along the spiral array of auditory nerve fibers. To achieve this, the processing channels of the sound processor; that is, specific frequency bands with their associated signal processing paths, are mapped to a set of one or more electrodes to stimulate a desired nerve fiber or nerve region of the cochlear. Such sets of one or more electrodes are referred to herein as “electrode channels” or, more simply, “channels.”

Conventional cochlear implants have limitations that may produce undesirable effects for recipients. One fundamental problem that limits the spatial resolution of multi-channel cochlear implants is referred to as “current spread” and is illustrated in FIG. 1. Although stimulation through one channel is intended to excite a single nerve region, in fact the actual locus of neural excitation can be broad and complex due the spread of current throughout the conducting fluids and tissues of the cochlea.

FIG. 1 is a graph illustrating the voltage created at various electrodes in response to the application of current at one electrode. Voltage profile 101 illustrates the voltage created at different nerve regions of the cochlear (“tissue voltage”) at a plurality of locations adjacent electrodes 104 of an electrode array 106 in response to current delivered to electrode number 11. Superimposed on FIG. 1 is a illustration of the current spread 102 emanating from the nerve region adjacent electrode 11 which causes voltage profile 101.

As illustrated by voltage profile 101, current delivered by electrode 11 may spread over a potentially wide spatial extent of neighboring nerve regions. This current spread may extend, for example, to nerve regions adjacent to distant electrodes 1 and 22 of the 22 electrodes of electrode array 106. As a result, a stimulating voltage 101 arises not only in the nerve region adjacent electrode 11 but also at more distant nerve regions in the tissue. As shown in FIG. 1, the stimulating voltage 101 is strongest or most intense near electrode 11, dropping off slowly and, in this example, remaining non-negligible at all regions in the cochlea nerve adjacent to electrode array 106. As a result, in addition to the nerve fibers adjacent electrode 11, other nerve fibers in the cochlea are stimulated by current applied to electrode 11. This may produce a distributed place-pitch perception, rather than the single pitch percept that was intended by the stimulating method.

This problem is exacerbated when current flows concurrently from two or more electrodes, as would occur when representing a sound with multiple frequency components. When two or more channels are activated concurrently, the locus of excitation is not the simple union of their individual loci because of the nonlinearity of the neural excitation process. Instead, neurons that fall outside of the individual loci (i.e. those which would not respond to any one channel) may nevertheless be excited by the summed current fields. This results in the well-known phenomenon of “channel interaction” or “channel overlap.” Channel interaction can lead to unpredictable loudness fluctuations, and smearing of the spatial representation of spectrum.

FIG. 2 shows the consequences of current spread 202 ₉ and 202 ₁₃ when current concurrently flows from two electrodes 9 and 13 of electrode array 106. Voltage profile 208 is generated in response to stimulating electrode 9, while voltage profile 210 is generated in response to stimulating electrode 13. Voltage profile 212 is the sum of voltage profiles 208 and 210; that is, voltage profile 212 is generated in response to simultaneously stimulating electrodes 9 and 13. As shown in FIG. 2, the combined currents produce a stimulus voltage in the nerve region adjacent to each electrode which is greater than intended, as well as a high voltage 214 in the nerve region between the electrodes 9 and 13.

This summation of stimulus voltages has many undesirable perceptual consequences, particularly when many electrodes are activated simultaneously to represent a complex sound with multiple frequency components. For example, such stimulation may result in unpredictable and excessive loudness and loss of spectral shape, that is, the peaks of the frequency-place profile are distorted by the summation of fields.

Almost all successful stimulation strategies in clinical use today circumvent channel interaction by using sequential pulsatile stimulation. Such strategies deliver stimulation through only one channel at any given instant. Stimulation is time-multiplexed across channels at rates high enough to produce a fused percept for the recipient. Although monopoles excite broad spatial extents of the nerve array, spatial/spectral information is nevertheless adequately conveyed, presumably by the trajectory of the centroids of those ranges. In this way a reasonable representation of a sound's time-varying magnitude spectrum can be appreciated by the recipient, such that formant peaks can be perceived.

Prior to the widespread adoption of sequential-monopolar stimulation for clinical use, several more complex channel configurations were explored with the objective of producing more focused electrical fields and hence narrower stimulation regions. These included bipolar stimulation with longitudinally- radially-, and obliquely-oriented dipoles, tripoles or quadrupoles, common ground, and complex multipolar channels derived by so-called “current deconvolution.” These more complex channel configurations have not resulted in speech understanding gains.

SUMMARY

In one aspect of the invention, a cochlear implant is disclosed, the cochlear implant configured to provide stimulation of a discrete stimulation region of a particular recipient's cochlear, said discrete stimulation region being defined by interference of stimulating and limiting signals simultaneously applied to electrode channels of an electrode array implanted in the cochlear, wherein the stimulating and limiting signals are determined based upon transimpedance measurements representing specific spread functions of the particular recipient.

In another aspect of the invention, a method for providing stimulation of a discrete stimulation region of a particular recipient's cochlear is disclosed. The method comprises: measuring recipient-specific current spread functions in the form of a matrix of transimpedance values between stimulated and idle electrodes; computing an inverse matrix of transadmittance values is computed using the measured transimpedance values; and utilizing the matrix of transadmittance values to determine the required vector of electrode currents that will produce the desired vector of stimulating voltages.

BRIEF DESCRIPTION OF THE DRAWINGS

Embodiments of the present invention will be described in conjunction with the accompanying drawings, in which:

FIG. 1 shows a graph illustrating the voltage generated in various regions in a tissue in response to the application of current to one tissue region;

FIG. 2 shows a “current spread” graph when two electrodes are stimulated and the corresponding electrode array;

FIG. 3 is a diagram of the components of an exemplary cochlear implant in which embodiments of the present invention may be advantageously implemented;

FIG. 4 is a flowchart of the operations performed in accordance with one embodiment of the present invention;

FIG. 5A shows a graph of the components of stimulating signal in accordance with an embodiment of the present invention;

FIG. 5B shows equations used to construct stimulating and limiting signals in accordance with an embodiment of the present invention;

FIG. 6A is a graph showing voltage profiles generating by applying the same current to each electrode of an electrode array;

FIG. 6B is a graph showing current levels required to apply a stimulating signal on one electrode while applying a limiting signal on the remaining electrodes in the electrode array; and

FIG. 7 is a graph showing the voltage profiles and current levels of a phased-array stimulation in accordance with an embodiment of the present invention.

DETAILED DESCRIPTION

Embodiments of the present invention are directed to a medical stimulation device such as a cochlear implant that, when operably implanted in a recipient's cochlear, provides stimulation of one or more spatially-restricted contiguous portion(s) of the spiral array of auditory nerve fibers in the cochlear, referred to herein as discrete stimulation regions. Each discrete stimulation region is defined by the constructive and/or destructive interference of stimulating and limiting signals simultaneously applied to electrode channels of an implanted electrode array, the stimulating and limiting signals being determined based upon transimpedance measurements of intracochlear electrode channels of the implanted electrode array representing specific spread functions of an individual recipient.

In certain embodiments the stimulating signal is applied through a targeted electrode channel; that is, one or more successive electrodes which is/are adjacent to the discrete stimulation region. The targeted electrode channel is selected to represent a sound based on the outputs of a sound processor to stimulate neural activity in the discrete stimulation region to thereby cause a percept of the represented sound. The size of the discrete stimulation region is defined by the limiting signal(s) applied to electrode channel(s) other than the targeted electrode channel, and which negate(s) current spread which would otherwise occur in response to the stimulating signal. Thus, electrodes of the electrode array may either apply a stimulating signal, no signal or a limiting signal.

Once having read the present application, one of ordinary skill in the art would appreciate that various combinations of stimulating and limiting signals may be used to negate the effects of current spread and to provide precise stimulation for spatial-temporal applications.

There are a variety of stimulating medical devices in which embodiments of the present invention may be advantageously implemented. FIG. 3 is a perspective view of an exemplary stimulating prosthetic hearing implant, a cochlear implant 300 in which an embodiment the present invention is implemented. The relevant components of outer ear 301, middle ear 305 and inner ear 307 are described next below, followed by a description of cochlear implant 300.

An acoustic pressure or sound wave 303 is collected by outer ear 301 (that is, the auricle) and channeled into and through ear canal 302. Disposed across the distal end of ear canal 302 is a tympanic membrane 304 which vibrates in response to acoustic wave 303. This vibration is coupled to oval window or fenestra ovalis 315 through three bones of middle ear 305, collectively referred to as the ossicles 317 and comprising the malleus 313, the incus 309 and the stapes 311. Bones 313, 309 and 311 of middle ear 305 serve to filter and amplify acoustic wave 303, causing oval window 315 to articulate, or vibrate. Such vibration sets up waves of fluid motion within cochlea 332. Such fluid motion, in turn, activates tiny hair cells (not shown) that line the inside of cochlea 332. Activation of the hair cells causes appropriate nerve impulses to be transferred through the spiral ganglion cells (not shown) and auditory nerve 338 to the brain (not shown), where they are perceived as sound.

Cochlear prosthesis 300 comprises external component assembly 342 which is directly or indirectly attached to the body of the recipient, and an internal component assembly 344 which is temporarily or permanently implanted in the recipient.

External assembly 342 typically comprises a sound transducer 320 for detecting sound, and for generating an electrical audio signal, typically an analog audio signal. In this illustrative embodiment, sound transducer 320 is a microphone. In alternative embodiments, sound transducer 320 may comprise, for example, more than one microphone, one or more telecoil induction pickup coils or other device now or later developed that may detect sound and generate electrical signals representative of such sound.

External assembly 342 also comprises a speech processing unit 316, a power source (not shown), and an external transmitter unit 306. External transmitter unit 306 comprises an external coil 308 and, preferably, a magnet (not shown) secured directly or indirectly to the external coil 308.

Speech processing unit 316 processes the output of microphone 320 that is positioned, in the depicted embodiment, by outer ear 301 of the recipient. Speech processing unit 316 generates coded signals, referred to herein as a stimulation data signals, which are provided to external transmitter unit 306 via a cable (not shown). Speech processing unit 316 is, in this illustration, constructed and arranged so that it can fit behind outer ear 301. Alternative versions may be worn on the body or it may be possible to provide a fully implantable system which incorporates the speech processor and/or microphone into the internal component assembly 344.

Internal components 344 comprise an internal receiver unit 312, a stimulator unit 326 and an electrode assembly 318. Internal receiver unit 312 comprises an internal transcutaneous transfer coil (not shown), and preferably, a magnet (also not shown) fixed relative to the internal coil. Internal receiver unit 312 and stimulator unit 326 are hermetically sealed within a biocompatible housing. The internal coil receives power and data from external coil 308, as noted above. A cable or lead of electrode assembly 318 extends from stimulator unit 326 to cochlea 332 and terminates in an array 334 of electrodes 336. Signals generated by stimulator unit 326 are applied by electrodes 336 to cochlear 332, thereby stimulating the auditory nerve 338.

In one embodiment, external coil 308 transmits electrical signals to the internal coil via a radio frequency (RF) link. The internal coil is typically a wire antenna coil comprised of at least one and preferably multiple turns of electrically insulated single-strand or multi-strand platinum or gold wire. The electrical insulation of the internal coil is provided by a flexible silicone molding (not shown). In use, internal receiver unit 312 may be positioned in a recess of the temporal bone adjacent to outer ear 301 of the recipient.

It should be appreciated that, as noted elsewhere herein, embodiments of the present invention may be implemented in stimulating prosthetic hearing implants other than cochlear implant 300. For example, while cochlear implant 300 is described as having external components, in alternative embodiments, cochlear implant 300 may be a totally implantable prosthesis. In one exemplary implementation, for example, sound processing unit 316, including microphone 320, a sound processor and/or a power supply may be implemented as one or more implantable components.

As shown in FIG. 3, cochlear implant 300 is further configured to interoperate with a wireless user interface 346 to facilitate implant configuration and control by the recipient, and an external processor 342 such as a personal computer, workstation or the like, implementing a hearing implant fitting system.

As one of ordinary skill in the art will appreciate, the present invention may be used in connection with any speech strategy now or later developed, including but not limited to, Continuous Interleaved Sampling (CIS), Spectral PEAK Extraction (SPEAK), and Advanced Combination Encoders (ACE™). An example of such speech strategies is described in U.S. Pat. No. 5,271,397, the entire contents and disclosures of which is hereby incorporated by reference herein. The present invention may also be used with other speech coding strategies now or later developed. Among other things, these strategies offer a trade-off between temporal and spectral resolution of the coded audio signal by changing the number of frequency channels chosen in the signal path.

As noted, channel interaction may limit the effectiveness of a cochlear implant. In particular, because of channel interaction, speech strategies which simultaneously stimulate multiple electrodes have only limited clinical success. Various methods have been devised over the years to address the problem of current spread to minimize the resulting channel interaction.

The most effective approach has been the adoption of so-called “sequential” stimulation strategies, wherein only one electrode produces current at any given moment. Since no two electrodes generate stimulating voltages simultaneously, the summation of FIG. 2 is not likely to occur. By rapidly sequencing the stimulating current among selected electrodes, a percept may be generated that is “apparently” simultaneous.

Sequential stimulation avoids the difficulties of loudness summation, but does not avoid the fundamental problem that a stimulus intended to stimulate only one pitch region, i.e. nerve region, may broadly stimulate many pitch regions. Thus, sequential stimulation may not produce excitation only at the desired region(s). Sequential stimulation also imposes serious timing constraints that prevent conveyance of some temporal structures contained in a sound. If two components of a sound would call for stimulation of two electrodes at the same moment in time, the sequential strategies forbid this. Harmonic relationships and other temporally-encoded cues may be seriously compromised as well. Although, sequential pulsatile strategies have been adopted to circumvent some of the worst effects of channel interaction, those strategies, along with the current spread itself, impose significant constraints on the temporal and spatial patterns of nerve activity that can be artificially produced with a cochlear implant.

Previous stimulating methods to overcome these limitations have not met with success. These methods are described in U.S. Pat. No. 4,648,403; Van Compernolle, Dirk, “Speech Processing Strategies for a Multichannel Cochlear Prosthesis,” Ph.D. Dissertation, Stanford University (1985); Van Compernolle, Dirk, “A computational model of the cochlea used with cochlear prosthesis patients,” in: Acoustics, Speech, and Signal Processing, IEEE International Conference on ICASSP '85, Volume 10, pp. 427-429 (1985); White, RL and Van Compernolle, D, “Current spreading and speech-processing strategies for cochlear prostheses,” Ann. Otol. Rhino. Laryng. 96 (Suppl. 128), 22-24 (1987); Townshend B, et al., “Pitch perception by cochlear implant subjects,” J. Acoust. Soc. Am., 82(1):106-115 (1987); Townshend B, and White R L, “Reduction of Electrical Interaction in Auditory Prostheses,” IEE Tran. Biomd. Eng. BME-34:891-897 (1987), the entire contents and disclosures of which are hereby incorporated by reference herein. These conventional stimulating methods are limited by psychophysical measures of current spread, which are tedious and time consuming, and may produce computed solutions requiring enormous current. A further significant drawback is that psychophysical measures of current spread provides only an imprecise measure of the physical parameters of interest.

Bipolar stimulation, where current flows from one electrode to another nearby electrode, partially reduces channel interaction, but at the expense of significantly increased current requirements. The resulting improvement in spatial sharpening has generally been modest.

Prior tripolar or quadrupolar electrode configurations have been described as a means for narrowing the stimulus area of an electrode. Certain embodiments of such configurations are described in Jolly C N, Spelman F A, Clopton B M, “Quadrupolar stimulation for Cochlear prostheses: modeling and experimental data,” IEEE Trans. Biomed. Eng. 43(8):857-865 (1996); Clopton & Spelman, “Electrode configuration and spread of neural excitation: compartmental models of spiral ganglion cells,” Ann. Otol. Rhinol. Laryngol. 166:115-118 (Suppl. 1995); Miyoshi, et al., “Proposal of a new auditory nerve stimulation method for cochlear prosthesis,” Artif. Organs 20:941-946 (1996); Kral, et al., “Spatial resolution of cochlear implants: the electrical field and excitation of auditory afferents,” Hear Res. 121:11-28 (1998); Townshend, et al., “Pitch perception by cochlear implant subjects,” J. Acoust. Soc. Am. 82(1):106-115 (1987), the entire contents and disclosures of which are hereby incorporated by reference herein. The above and other prior art tripolar/quadrupolar approaches utilize fixed weights based upon mathematical models, or physiological measurements, and do not contemplate consideration of the individual recipients. Further, conventional tripolar/quadrupolar approaches generally ignore contributions from other more distant electrodes thereby limiting their ability to sharpen the field.

Embodiments of the present invention each overcome one or more of the above and other constraints and may allow any spatio-temporal pattern of excitation, within the spatial resolution of an electrode array. Embodiments of the present invention stimulate discrete nerve regions through the constructive and destructive interference of stimulating signals and interfering limiting signals, which may be of the same or opposite polarity, applied via multiple electrode channels of the electrode array. This interference generates a stimulating voltage at a discrete stimulation region within the cochlea. The stimulating and limiting signals are determined based upon transimpedance measurements from electrodes within the cochlea.

One advantage of this stimulation method is that the implemented strategy is customized to the recipient to improve his/her hearing response. Also, embodiments of the present invention may provide a substantial improvement in sound representation, better speech understanding, and improved recognition of musical information, such as melody and harmonic relationships.

In one embodiment of the present invention every electrode in the electrode array may apply a stimulating or limiting signal. Some electrodes might be used to select the current balance that minimizes spread of current with a limiting signal. Through the use of voltages at every electrode, the discrete nerve region may be precisely and simultaneously controlled using proper sources and sinks at each electrode. Thus, a recipient may receive a customized solution for determining the necessary current weights to improve hearing.

In certain embodiments, the stimulation involves a contribution from every electrode to every voltage. Further embodiments may apply a similar concept to a subset of electrodes. For certain applications and implants, such further embodiments may be advantageous, because fewer current sources and/or arithmetic operations may be required. For example, tripolar stimulation has previously been suggested as a possible means of ameliorating current spread. In a tripolar configuration, some of the current flowing from a center electrode returns to each of two neighboring electrodes, and some may return to a distant electrode, i.e. an electrode outside of the cochlea. Tripolar configurations described in the literature involve fixed weights, such as 30% of the current to each neighbor, and 40% to the distant electrode. But due to anatomical and positioning variations the optimal weights for minimizing current spread with a tripole would generally be different for each group of three electrodes. In contrast, embodiments of the present invention having phased-array stimulation are not subject to such restrictions, and may be used to customize the individual the weights for each tripolar group of electrodes. Other embodiments combine positive and negative currents of appropriate magnitude through more than one electrode to create a desired profile of stimulating voltages across N electrode sites.

FIG. 4 is a flow chart of the operations performed in accordance with one embodiment of the present invention. The operations are described next below with reference to additional FIGS. 5A and 5B. At block 402, recipient-specific current spread functions are measured in the form of a matrix of transimpedance values between stimulated and idle electrodes.

For each intracochlear electrode, a monopolar biphasic pulse is delivered in one embodiment of the present invention. When a monopolar current is passed through an intracochlear electrode 504 a corresponding voltage can be measured at any other intracochlear electrode 504. The transimpedance between the two electrodes is defined as the ratio of the measured voltage to the delivered current. For frequencies of interest in cochlear applications, the cochlear fluids and tissue are essentially resistive such that the voltage is approximately instantaneously proportional to the current to a close approximation. Thus the reactive component of the transimpedance may be ignored. Nevertheless, the term transimpedance is used throughout this application in preference to transresistance as a reminder that the analysis presented herein may be carried out with complex arithmetic in those applications in which reactive components are nonnegligible.

For each intracochlear electrode 504, a monopolar biphasic pulse is delivered at the maximum comfortable current level. Return current is measured, for example, at one or more extra-cochlear electrodes, the arrangement of which is well-known in the art. The voltage pulse at each of the remaining electrodes 504 is measured. In one embodiment, the phase width is selected to be long enough for the voltage pulse waveform to plateau, but short enough to permit relatively high currents while remaining comfortable, thus maximizing the signal-to-noise ratio of the voltage measurement. Based on the applied current and measured voltage, the transimpedance of each combination of electrodes is determined as described below.

Referring to FIG. 5A, each current I₁ through I₂₂ flowing simultaneously through each of the 22 different electrodes 504 results in a current spread 510 ₁ through 510 ₂₂, respectively. The instantaneous voltage 512 V₁ . . . V₂₂ adjacent to each electrode 504 is the sum of 22 separate components 514 ₁ through 514 ₂₂. Each component 514 is proportional in magnitude and sign to the current from one of the electrodes 504 I₁ . . . I₂₂ (in units of mA), represented by the corresponding measured transimpedance value, z₁ . . . Z₂₂ (in units of V/mA (kΩ)). Thus the stimulating voltage 512 ₁₁ at the site of electrode 11 may be expressed as a weighted sum of the weighted currents through all stimulating electrodes, as shown in FIG. 5A and Equation (1) below: V ₁₁ =Z ₁ I ₁ +Z ₂ I ₂ + . . . +Z ₂₂ I ₂₂  (4)

A similar equation may be written for the voltage adjacent to each of the other electrodes 504 as shown in FIG. 5B. Thus, in a 22 electrode system, 22 simultaneous equations may be written which describe the voltage applied to the nerve region adjacent to each electrode as a weighted sum of the same set of 22 currents I₁ . . . I₂₂. The set of 22 simultaneous equations in FIG. 5B, where each weight, transimpedance Z, now has two subscripts, one indicating the associated current, and one indicating the electrode site whose voltage is being summed. The set of simultaneous equations in FIG. 5B, may be represented in vector/matrix notation as shown in Equation (2): v=Zi_(e)  (2) where I_(e) represents the column vector of 22 currents I₁ . . . I₂₂ flowing through the electrode, Z represents the square matrix of weights z_(1,1) . . . z₂₂, and V represents the column vector of 22 stimulating voltages V₁ . . . V₂₂ at the discrete stimulation regions adjacent the corresponding electrodes I₁ . . . I₂₂.

As noted, both the currents and voltages are represented as functions of time, indicating that the matrix equation represents an instantaneous calculation. It should be appreciated, however, that as noted above, the preferred embodiment uses simple scalar values. This reflects the assumption that voltages are instantaneously proportional to currents (equivalent to assuming that the tissue impedances is purely resistive with no reactive component). So if a current changes, all resulting voltages change instantaneously and proportionally. While this is a relatively accurate approximation, it is not exactly correct. At the expense of greatly increased computational burden, embodiments of the invention may be practiced using complex values in the matrices, recognizing the reactive (non-resistive) part of the relationship between currents and voltages. Then instead of being functions of time, the variables in equations 1 and 2 would be functions of the Laplace transform variable s.

For a stimulating electrode j and measuring electrode i, transimpedance z_(ij) has units of V/mA (kΩ) and is given by Equation (3): z _(ij) =v _(i) /i _(j)  (3) where, v_(i) is the measured peak voltage at measuring electrode i, and i_(j) is the amplitude of the current pulse applied to electrode j. In this way a 22×22 matrix Z_(m) of transimpedance values may be determined for each recipient, as shown in Equation (4):

$\begin{matrix} {Z_{m} = \begin{bmatrix} z_{1,1} & z_{1,2} & \ldots & z_{1,22} \\ z_{2,1} & z_{2,2} & \ldots & z_{2,22} \\ \vdots & \vdots & \ddots & \vdots \\ z_{22,1} & z_{22,2} & \ldots & z_{22,22} \end{bmatrix}} & (4) \end{matrix}$

Each column p of the matrix represents the spread function for stimulation through electrode p, with a peak at the diagonal value z_(pp). The values on the diagonal of Z_(m) typically can not be measured explicitly due to polarization effects. A current-carrying electrode is polarized by electrochemical gradients that arise across the metal/electrolyte boundary. This precludes using the same electrode both to deliver current and to measure potential in the cochlear fluid. Instead, values on the diagonal may be extrapolated, preferably using the highest slope among adjacent pairs in the same row and column to avoid underestimating the sharpness of the spread function, as this would result in unnecessarily high peak currents in the focused stimuli. The effect of errors among the extrapolated values is considered below.

As one of ordinary skill in the art would appreciate, the matrix Z_(m) is very nearly diagonally symmetric. For example, for a three-port network where a single node serves as both the return current path and the voltage measurement reference, the reciprocity theorem holds that diagonal symmetry obtains—i.e. z_(ij)=z_(ji) for all i and j.

In fact, due to the noted polarization phenomenon, the extracochlear electrodes could not serve as both return current path and voltage reference. Instead, a surface electrode on the recipient may serve as the reference. However the observed deviations from diagonal symmetry were small and comparable to the noise in the measurements. This implies that the tissue adjacent the return electrode is essentially neutral, or equipotential with the external electrode, and that the voltage gradients of consequence occur within and about the cochlea. Therefore it may be assumed that deviations from diagonal symmetry are dominated by noise in the measurement. In order to reduce this noise the diagonally opposite elements may be averaged to compute a new transimpedance matrix, as shown in Equation (5):

$\begin{matrix} {Z = {\frac{1}{2}\left( {Z_{m} + Z_{m}^{T}} \right)}} & (5) \end{matrix}$ where Z_(m) ^(T) denotes the transpose of Z_(m) (matrix Z is diagonally symmetric by definition).

Certain embodiments take advantage of the Reciprocity Theorem to avoid measuring the entire forward matrix [z], omitting measurement of cells above the main diagonal of matrix [z]. The unmeasured values are then filled in by transposing values from below the main diagonal. Alternatively, embodiments may omit the measurement of cells below the main diagonal, and fill these measurements by transposing values from above. In either case such embodiments may reduce the number of measurements needed to determine matrix [z] by half.

FIG. 6A is a graph showing a set of voltage profile curves 602 ₁ . . . 602 ₂₂ arising from stimulation of one electrode 604 at a time from an electrode array 606. The vertical values, in V/mA, of the points on each curve 602 represent one column of the transimpedance matrix [z] shown in FIG. 5B. Curve 602 ₁ represents the first column of FIG. 5B, curve 602 ₂ the second column, curve 602 ₃ the third column, and so on. Each voltage profile curve 602 is produced by passing 1 mA through each electrode 604. Curves 602 may have differing peak values and widths depending on a variety of factors including the type of electrode array used, how the electrode array is implanted and the recipient.

The forward matrix [z] may also be determined as follows. A fixed amplitude current waveform is applied to a single electrode (number j) with respect to a distant electrode. The corresponding voltage waveform may be measured at each of the other unstimulated electrodes with respect to a non-current carrying reference electrode at a neutral site in the tissue. Signal averaging may be used as necessary to obtain a clean measurement. Some feature of the stimulus waveform is selected for computing the ratio of voltage to current. For example, if the stimulus current is a rectangular biphasic pulse, the peak of the first phase could serve as a reference feature. The value for z_(ij) would be equal to the ratio of the peak voltage on electrode i to the peak current flow through electrode j. This set of measures would be used to define one column of [z] except for the entry z_(jj), the weight for the stimulated electrode, for the reasons noted above. The missing value of z_(jj) may be estimated by extrapolating from curve(s) fitted through several values above and/or below electrode j. The entire process above may be repeated for each stimulating electrode to populate the entire forward matrix [z].

The required voltage measurements may be carried out in any location internal or external to cochlear implant 100, such as in sound processor 316, computer 342, etc., and the transimpedance values or representations thereof may be transmitted via to sound processor 316 and stimulator unit 326 as necessary.

Returning to FIG. 4, at block 404, an inverse matrix of transadmittance values is computed using the computed transimpedance values. The matrix of transadmittance values is utilized to determine the required vector of electrode currents that will produce the desired vector of stimulating voltages. As will be described below, each column of the transadmittance matrix comprises a set of numerical weights (transadmittance values) defining the current contribution from each electrode that is required to produce a non-zero intrascalar voltage at a single discrete stimulation region. As such, each such vector of weights defines an electrode channel.

Because Z is square, and provided that it is nonsingular, we can compute its inverse Y=Z⁻¹. Then solving equation (3) for i_(e) yields equation (6): i_(e=Z) ⁻¹v_(d)=Yv_(d)  (6) The elements of Y are transadmittance values with units of mA/V (millimhos).

Equation 6 provides a solution to the inverse problem and from this equation, for any instantaneous set of desired stimulating voltages [V] the instantaneous currents [I] necessary to generate the stimulating voltages may be computed. Embodiments of the present invention may compute that the inverse matrix [z] once, after which it may be used continuously to compute any desired series of currents [I].

The constraint that Z be non-singular is not mathematically problematic, as a singular matrix will not arise with actual weighting constants for an actual electrode set. In practical terms, the corresponding problem is that, if the profiles defined by the rows of Z are quite broad and almost horizontal, its Condition Number may be large, such that the required currents may also be very large.

The graph shown in FIG. 6B shows part of the inverse matrix Y from Equation 6. Each curve 610 represents one column of matrix Y. The vertical values have units of mA/V. FIG. 6B shows the current required from each electrode 604 on an electrode array 606 to generate a stimulating voltage of 1.0 V at a single electrode site, and 0 V at all other electrode sites. For clarity, FIG. 6B shows five of the possible 22 curves. For example, curve 610 ₇ illustrates the currents required on each electrode to generate a stimulating voltage of 1.0 V adjacent to electrode 7 (and 0 V at all other sites). Similarly, curve 61 ₁₁ shows the currents required to generate a 1.0 V stimulating voltage only adjacent to electrode 11.

Producing a specified voltage v_(d) _(p) at a single discrete stimulation region, p, and exactly zero volts at every other region along the cochlear represents optimally focused single-place stimulation. The required vector of currents is the product of scalar v_(d) _(p) and column p of Y:

$\begin{matrix} {i_{e} = {v_{d_{p}}\begin{bmatrix} y_{1p} \\ \vdots \\ y_{22p} \end{bmatrix}}} & (7) \end{matrix}$

The transadmittance elements (y_(1p) . . . y_(22p)) of column p constitute a set of 22 signed weights that define the current ratios needed to produce a stimulating voltage only at place p. Thus each column of Y defines a multipolar electrode channel; that is, a set of electrodes among which currents are passed with specified ratios, in order to excite nerve fibers at a single spatially-restricted discrete stimulation region within the cochlea.

Embodiments of the present invention may generate a set of voltages across a 22 electrode array by computing the currents necessary to generate each voltage individually, and then summing the required currents on each electrode. This is equivalent to computing the current vector [I] in Equation 6.

FIG. 7 shows two graphs that illustrate the voltage profile curves 702 and current levels 708 for each electrode 704 on an electrode array 706 to produce 1.0 volt at electrodes 9 and 13, and zero volts on the remaining electrodes 704. On electrode 9 there is a positive current level 708 ₉ corresponding to a stimulating signal. While on the adjacent electrodes 8 and 10 there is a negative current 708 ₈ and 708 ₁₀ corresponding to a limiting signal. On electrode 20 there is a slight positive current 708 ₂₀ which has a limiting signal providing interference.

Net voltage curve 730 is the sum of all of voltage profile curves 702. Net voltage curve 730 represents the stimulating voltage profiles 702 that would arise from the combination of current levels 708. Various combinations of current levels may be generated by embodiments of the present invention to produce desired net voltage to stimulate discrete nerve regions without creating channel interference. Simultaneous stimulation may be achieved by such embodiments.

FIG. 7 also shows the advantages the present invention in preventing channel interference when compared with FIG. 2. In both cases the strategy sought to stimulate only two regions. In the simple case of FIG. 2 the voltage profile was quite different from the intended stimulation. Due to current spread, stimulating voltages were created across the entire length of the electrode array. However by applying the phased-array currents of the present invention, shown in FIG. 7, the profile of stimulating voltages may be precisely achieved.

As is illustrated in FIG. 7, the peak currents required in a phased-array stimulus may be substantially higher than those used for monopolar stimulation. In practice, the broader the single-electrode profiles, shown in FIG. 6A, the higher will be the required currents. As the single-electrode profiles broaden, the Condition Number of the matrix [z] increases, resulting in higher peak current requirements. Higher currents also imply greater power dissipation within the implanted package. Also, if the curves are very flat and the Condition Number of matrix [z] is large, the method demands greater precision in the arithmetic and current sources, and becomes sensitive to small drifts in the coefficients of [z] and [Y]. For all of these reasons, embodiments of the present invention using phased-array stimulation may be implemented in Nucleus Contour™ or similar implant devices having reasonably sharp single-electrode profiles. Also, when electrodes are positioned along the lateral wall of the cochlea larger current, power and precision requirements may be required.

The values of coefficients in the forward matrix [z] may vary slowly over time due to scar tissue formation around the electrode, bone growth, electrolyte shifts, etc. In that case it may be necessary to remeasure and update the matrices from time to time, such as when visiting a clinic or doctor. It might be possible, for example, for the implant system to remeasure [z] each time the system is turned on.

With existing implant designs, if a current source saturates and is unable to deliver the intended current, the result is benign in that it produces a sound that is less loud than intended. In embodiments of the present invention, the phased-array stimulation may produce loud and possibly uncomfortable sounds when the current source saturates. This may occur when the saturated source fails to cancel in part the current from another unsaturated source. Such embodiments may take precautions by designing an implant that prevents this from occurring, or scales all current sources down when any electrode approaches saturation.

Implementation of a phased-array stimulus of embodiments of the present invention may depend on the presumption that the fluid in the cochlea and tissue are electrically linear and approximately resistive within the frequency range of interest, i.e. up to the highest meaningful spectral component of the stimulus current waveforms. Current preliminary animal and human data support this presumption. As one of ordinary skill in the art would appreciate, alternative embodiments of the present invention may be modified to compensate for non-linear responses of the tissue to changes in the fluid in the cochlea.

In the above embodiments, the stimulation involved a contribution from every electrode to every voltage. A phased-array stimulation of the present invention may also be applied to any subset of M electrodes out of the total N electrodes in an array. For example, with a 22-electrode array, 16 different 7-electrode groups could be formed at one-electrode intervals along the array. The electrodes in any subset need not be contiguous or uniformly spaced. Various combinations of subsets of electrodes are envisaged by the present invention.

Implementations of the present invention may use an imperfect or suboptimal solution when generating the stimulating and limiting signal. An “imperfect solution” or “suboptimal solution,” may refer to phase array stimulation that sets the current within a fixed level or constraints. Such a solution may be completely satisfactory for certain recipients and applications, and may provide some practical economy in the implementation. There are two situations where an imperfect solution may be desired. First, if some elements of the matrix [z] are sufficiently small or negligible, these elements may be arbitrarily set to zero without substantially influencing the resulting voltage profile. Such embodiments may reduce computational burden or the number of current sources required. Second, there may be instances where a perfect solution is impractical because such a solution requires peak currents in excess of the implant's capabilities or safety boundaries. In these cases, a suboptimal solution may be created within some set of constraints. For example, the largest current may be less than or equal to 2 mA, and the maximum number of non-zero weights in a column may be less than 10. The suboptimal solution matrix [z] could be designed by optimizing some metric of merit, e.g. mean square deviation from target voltage profile, while meeting the specified constraint.

The present invention may be implemented in implant devices where there are more channels than electrodes. The discussion above describes methods for computing a perfect solution or imperfect solution for N stimulation voltages using N electrodes. Embodiments of the present invention may also generate from the original forward measures [z] an N×L pseudoinverse matrix which would map N currents to L stimulation sites, where L is greater than N. This pseudoinverse may be generated by a least squares fit or other interpolation process beginning with the original trans-impedance measures of [z]. A pseudoinverse is described by Van Compernolle, Dirk, “A Computational Model of the Cochlea used With Cochlear Prosthesis Patients,” in: Acoustics, Speech, and Signal Processing, IEEE International Conference on ICASSP '85, Volume 10, pp. 427-429 (1985), the entire contents and disclosures of which are hereby incorporated by reference herein. In addition, Rodenhiser and Spelman, “A Method for Determining the Driving Currents for Focused Stimulation in the Cochlea,” IEEE Trans. Biomed. Eng. 42(4) (April 1995) pg. 337-342, which are hereby incorporated by reference herein, described a method for creating a pseudoinverse matrix similar to that of Van Compernolle. However Rodenhiser and Spelman did not teach the use of voltage measurements from unstimulated electrodes to determine the matrix. Instead, Rodenhiser and Spelman proposed use of a mathematical model of fields within a generic cochlea. Such a model approach would not be useful in practice due to individual differences among electrode positions and anatomy of the cochlea across implant recipients.

Matrix Y defines channel weights for the 22 intracochlear electrodes, but each channel also employs extracochlear current as well. The total current flowing into the tissue must equal the total current flowing out. Therefore the extracochlear current must be equal and opposite to the algebraic sum of the intracochlear currents. Its weight y_(x) _(p) for channel p is given by Equation (8):

$\begin{matrix} {y_{x_{p}} = {\underset{i}{- \sum}y_{ip}}} & (8) \end{matrix}$

When comparing stimulation at two places k and p, the fact that |v_(k)|>|v_(p)| does not imply that the stimulus at k is functionally stronger than that at p because the magnitude of the biological response depends upon other factors as well (e.g. proximity of contact to the neurons, local neural density). In practice, stimulus intensity for each place is specified independently with respect to some place-specific behavioral criterion level such as detection threshold. We therefore normalized the weights for each channel to produce a dimensionless set of weights w_(p)=(w_(1p) . . . w_(22p)) as follows:

$\begin{matrix} {\begin{bmatrix} w_{1p} \\ \vdots \\ w_{22p} \end{bmatrix} = {\frac{1}{y_{pp}}\begin{bmatrix} y_{1p} \\ \vdots \\ y_{22p} \end{bmatrix}}} & (9) \end{matrix}$ The extracochlear current weight is similarly normalized as w_(x) _(p) =y_(x) _(p) /y_(pp).

As expected, for each channel p diagonal element y_(pp) was invariably the largest of the original channel weights. Thus the largest normalized weight w_(pp) was equal to 1.0. Normalization of the weights preserves the ratios of currents used for stimulation through a single channel p, so that it still produces exactly 0V at all other places k≠p.

The normalized channel weights form the columns of a new matrix W that was used to compute current vectors for experimental stimulation. The input was no longer a vector of desired voltages. Instead the currents were computed from Equation (10): i_(e=Wi) _(c)  (10) where input i_(c) was a vector of stimulus intensities expressed as “on-center” currents (in mA). The on-center current for a channel is its contribution to the current carried by the electrode located at its stimulation place, while it contributes smaller currents of equal or opposite polarity to the other “off-center” electrodes. Specifying a channel's stimulus intensity in mA this way has the advantage that it allows intuitive comparison with monopolar channels that deliver on-center current only.

Equation (9) provides an instantaneous solution to the inverse problem. In practice each of the stimulus intensities (i_(c) ₁ (t) . . . i_(c) ₂₂ (t)) is a function of time. For example those functions might be a set of biphasic pulse trains with varying intensities and pulse widths, or sinusoids of varying frequency and amplitude. Consequently the electrode currents (i_(e) ₁ (t) . . . i_(e) ₂₂ (t)) are also functions of time. A new current vector i_(e) must be computed each time that any signal i_(c) _(p) (t) changes value. In general for sampled signals, a new matrix multiplication is required for every sample interval. Thus by applying electrode current waveforms i_(e)(t) derived from repeated solution of equation (10) it is possible to deliver 22 completely independent stimulus waveforms to each of the 22 intracochlear places.

As noted the diagonal elements of matrix Z were necessarily estimated by extrapolation from surrounding values. It is important to understand the extent to which estimation errors compromise the resulting stimulation voltages.

If Z is the true transimpedance matrix, we define the imperfect matrix Z′ as Z′=Z+Δ  (11) where Δ is the diagonal matrix of estimation errors δ_(pp)=z′_(pp)−z_(pp). Using imperfect transimpedances Z′ and given a desired vector of stimulating voltages v_(d) the solution of equation (5) yields an imperfect vector i′_(e) of stimulating currents i′_(e)=Z′⁻¹v_(d)  (12)

The question at hand is this: “When the imperfect currents of vector i′_(e) are passed through the electrodes how will the actual voltages v′ differ from the desired voltages V_(d)?” From equation (2) the actual voltages are given by v′=Zi′_(e)  (13)

This may be compared with the correct voltages v_(d) that would have arisen if the transimpedances really were given by Z′. From equation (12), then: v_(d)=Z′i′_(e)  (14)

Subtracting equation (14) from equation (13) yields the vector of error voltages ε:

$\begin{matrix} \begin{matrix} {ɛ = {v^{\prime} - v_{d}}} \\ {= {{Zi}_{e}^{\prime} - {Z^{\prime}i_{e}^{\prime}}}} \\ {= {\left( {Z - Z^{\prime}} \right)i_{e}^{\prime}}} \\ {= {{- 1}\Delta\; i_{e}^{\prime}}} \end{matrix} & (15) \end{matrix}$

Two observations from equation (15) are noteworthy. First, since Δ is diagonal, the error voltage ε_(p) at place p is determined entirely by the corresponding estimation error δ_(pp) and the current through electrode p: ε_(p)=−δ_(pp) i′ _(ep)  (16)

Second, the negative multiplier of equation (16) indicates that overestimating the transimpedance z_(pp) results in a weaker than intended voltage and vice versa.

There are many advantages due to embodiments of the present invention. For example, a renewed interest in stimulus focusing has developed in recent years, motivated in part by the demands of strategies that seek to preserve within-channel fine temporal patterns of stimulation. These strategies deliver stimuli at instants in time determined from features of a processing channel's audio filter output (e.g., peaks or zero crossings). In general, those instants in time are asynchronous across channels, and so cannot satisfy the nonsimultaneity constraint of sequential stimulation. As a result, stimulus pulses inevitably will overlap in time resulting in overlapping stimulus regions which corrupt the temporal patterns. That is, due to the overlapping stimulus regions, any given neuron will generally be driven by pulses from multiple channels, receiving the composite of multiple temporal patterns rather than one. The narrowing of the stimulation regions minimizes or prevents this from occurring.

Another advantage of the present invention is the improved representation of periodicity (musical) pitch by an implementing cochlear implants. Besides its importance for music, pitch recognition is important for segregating one voice from multiple competitors and for transmission of semantic information in tonal languages. Perhaps ironically, periodicity pitch is only weakly conveyed by the repetition rate of an electrical stimulus, and then only for an isolated single frequency (as opposed to a chord) and over a limited frequency range. There is growing evidence that recognition of periodicity pitch by the central nervous system may depend upon detecting fine spatial patterning of harmonic complexes. The significantly improved focusing of the stimulation enables an implementing cochlear implant to recreate these detailed spatial patterns.

Although the present invention has been fully described in conjunction with several embodiments thereof with reference to the accompanying drawings, it is to be understood that various changes and modifications may be apparent to those skilled in the art. For example, in one embodiment, the limiting signals have a current which has a polarity opposite to that of the stimulating signal. In other embodiments, the limiting signal may also have a current of equal polarity such that a desired total net stimulation occurs at the discrete region. As another example, the embodiments of the present invention the electrodes are evenly spaced. In alternative embodiments, however, the electrodes have non-uniform spacing. Such changes and modifications are to be understood as included within the scope of the present invention as defined by the appended claims, unless they depart therefrom.

All documents, patents, journal articles and other materials cited in the present application are hereby incorporated by reference herein. 

1. A method for providing stimulation of one or more discrete stimulation regions of a particular recipient's cochlea via a cochlea implant, comprising: measuring recipient-specific current spread functions in the form of a matrix of transimpedance values of electrode channels of an electrode array implanted in the particular recipient's cochlea and having electrodes mounted thereon; calculating values for a diagonal of said matrix of transimpedance values using said measured recipient-specific current spread functions; computing a matrix of transadmittance values using said matrix of measured transimpedance values; determining, using said matrix of transadmittance values, a vector of electrode currents that will produce a desired vector of stimulating voltages resulting in stimulation of said one or more discrete stimulation regions; and applying stimulation, through at least one of said electrodes, to the recipient's cochlea using said vector of electrode currents.
 2. The method of claim 1, wherein measuring said matrix of transimpedance values comprises: simultaneously delivering a current through a desired number of said implanted electrodes; measuring the instantaneous voltage at each said electrode channel resulting from said current delivery, wherein said voltage at each electrode channel is a weighted sum of the currents through said desired number of electrode channels, and wherein each current is multiplied by a weight factor equivalent to a transimpedance between said electrode channel being measured and the particular electrode channel receiving the current being weighted; creating simultaneous voltage equations for each of said measured voltages; extracting said weight factors from said voltage equations; and inserting said extracted weight factors into a matrix, wherein each row of said matrix is an ordered collection of the weight factors for a particular voltage equation, and each column of said matrix represents the current spread function for stimulation through a particular electrode channel.
 3. The method of claim 2, wherein said desired number of electrodes comprise 22 electrodes, and wherein said matrix of transimpedance values (Z_(m)) is given as: $Z_{m} = \begin{bmatrix} z_{1,1} & z_{1,2} & \Lambda & z_{1,22} \\ z_{2,1} & z_{2,2} & \Lambda & z_{2,22} \\ M & M & O & M \\ z_{22,1} & z_{22,2} & \Lambda & z_{22,22} \end{bmatrix}$ wherein z_(1,1) through Z,_(22,22) are said weight factors.
 4. The method of claim 2, wherein said creating voltage equations comprises: using simplified computations that assume that the cochlea tissue of the particular recipient is purely resistive, thereby resulting in the use of non-perfect weight factors in said matrix of transimpedance.
 5. The method of claim 2, wherein creating said voltage equations comprises: using complex values recognizing the reactive (non-resistive) components of the cochlea tissue of the particular recipient, thereby resulting in approximately perfect weight factors in said matrix of transimpedance.
 6. The method of claim 2, wherein said matrix of transimpedance is measured for half of said matrix, and the other half of said matrix is calculated using the reciprocity theorem.
 7. The method of claim 1, wherein measuring said matrix of transimpedance values comprises: applying a fixed amplitude current waveform through a single electrode channel with respect to an extracochlear electrode of said electrode array; measuring the resulting voltage at each of a desired number of other electrode channels with respect to a non-current carrying reference electrode of said electrode array; calculating a ratio of each said measured voltage to said current applied through said single electrode channel; populating a column of a matrix with said calculated ratios; and repeating said above steps for each of said desired electrode channels.
 8. The method of claim 7, wherein calculating said ratio of each said measured voltage to said applied current comprises: using simplified computations that assume that the cochlea tissue of the particular recipient is purely resistive, thereby resulting in the use of non-perfect values in said matrix of transimpedance.
 9. The method of claim 7, wherein calculating said ratio of each said measured voltage to said applied current comprises; using complex values recognizing the reactive (non-resistive) components of the cochlea tissue of the particular recipient, thereby resulting in approximately perfect values in said matrix of transimpedance.
 10. The method of claim 7, wherein said desired number of electrodes comprise 22 electrodes, and wherein said matrix of transimpedance values (Z_(m)) is given as: $Z_{m} = \begin{bmatrix} z_{1,1} & z_{1,2} & \Lambda & z_{1,22} \\ z_{2,1} & z_{2,2} & \Lambda & z_{2,22} \\ M & M & O & M \\ z_{22,1} & z_{22,2} & \Lambda & z_{22,22} \end{bmatrix}$ wherein non-diagonal values z_(1,2) through Z_(22,21) are said calculated ratios of measured voltage to applied current for each electrode combination; and wherein the values for the diagonal, z_(i,i), of said matrix of impendance values (Zm) are calculated by extrapolating the values for the diagonal from a plurality of the calculated ratios .
 11. The method of claim 7, wherein said matrix of transimpedance is measured for half of said matrix, and the other half of said matrix is calculated using the reciprocity theorem.
 12. The method of claim 1, wherein computing said matrix of transadmittance values comprises: taking the inverse of said matrix of measured transimpedance values.
 13. The method of claim 1, wherein the transimpedance matrix is measured during programming by a clinician, and remeasurement is performed only upon visit to a clinician or doctor.
 14. The method of claim 1, wherein the transimpedance matrix is remeasured by the system each time the system is turned on.
 15. The method of claim 1, wherein there are more electrode channels than electrodes.
 16. The method of claim 1, wherein measuring said matrix of transimpedance and calculating said matrix of admittance are performed within said implant.
 17. A cochlear implant comprising: means for measuring recipient-specific current spread functions in the form of a matrix of transimpedance values of electrode channels of an electrode array implanted in the particular recipient's cochlea; means for calculating values for a diagonal of said matrix of transimpedance values using said measured recipient-specific current spread functions; means for computing a matrix of transadmittance values using said matrix of measured transimpedance values; and means for determining, using said matrix of transadmittance values a vector of electrode currents that will produce a desired vector of stimulating voltages resulting in stimulation of said one or more discrete stimulation regions.
 18. The implant of claim 17, wherein said means for measuring said matrix of transimpedance values comprises: means for simultaneously delivering a current to a desired number of said electrode channels; means for measuring the resulting instantaneous voltage at each said electrode channel, wherein said voltage at each electrode channel is a weighted sum of the currents through said desired number of electrode channels, and wherein each current is multiplied by a weight factor equivalent to the transimpedance between said electrode channel being measured and a particular electrode channel receiving the current being weighted; means for simultaneously creating voltage equations for each of said measured voltages; means for extracting said weight factors from said voltage equations; and means for inserting said extracted weight factors into a matrix, wherein each row of said matrix is an ordered collection of the weight factors for a particular voltage equation, and each column of said matrix represents the current spread function for stimulation through a particular electrode.
 19. The implant of claim 17, wherein means for measuring said matrix of transimpedance values comprises: means for applying a fixed amplitude current waveform to a single electrode channel with respect to an extracochlear electrode of said electrode array; means for measuring the resulting voltage at each of a desired number of other electrode channels with respect to a non-current carrying reference electrode of said electrode array; means for calculating a ratio of said each measured voltage to said current applied to said single electrode channel; means for populating a column of a matrix with said calculated ratios; and means for repeating said above steps for each of said desired electrode channels.
 20. The implant of claim 17, wherein means for computing said matrix of transadmittance values comprises: means for taking the inverse of said matrix of measured transimpedance values.
 21. The method of claim 1 wherein calculating values for a diagonal of said matrix of transimpedance values comprises: extrapolating the values for the diagonal from one or more curves fitted through a plurality of non-diagonal values of the transimpedance matrix.
 22. The implant of claim 17 wherein the means for calculating values for a diagonal of said matrix of transimpedance values comprises: means for extrapolating the values for the diagonal from one or more curves fitted through a plurality of non-diagonal values of the transimpedance matrix.
 23. The method of claim 1, wherein applying stimulation comprises: applying stimulation using a plurality of non-contiguous electrodes of the plurality of electrodes.
 24. The implant of claim 17, wherein the desired vector of stimulating voltages comprises stimulating voltages for application using a plurality of non-contiguous electrodes of the plurality of electrodes. 